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SUMMARY 


Closed  form  solutions  are  obtained  for  the  problem  of  a crack 
propagated  by  a 2-dimensional  rigid  wedge  of  finite  thickness  in  an 
infinite  power  law  viscoelastic  body.  The  distance  from  the  crack  tip 
to  the  point  of  contact  with  the  wedge  is  determined  from  the 
assumption  that  in  the  vicinity  of  the  crack  tip,  forces  of  cohesion 
act  to  produce  a smoothly  closing  crack  with  continuous  stresses. 
Finally,  a simple  formula  expressing  the  stress  intensity  factor  as 
a function  of  the  speed  of  penetration  and  the  length  of  the  free  crack 
surface  is  exhibited. 


INTRODUCTION 


In  a previous  paper  by  Walton,  Nachman  and  Schapery  [1],  hereafter 


i'  referred  to  as  WNS,  simple  closed  form  solutions  were  obtained  for  the 

i 2-dlmenslonal  problem  of  a rigid  indentor  sliding  over  a power  law 

!!  viscoelastic  halfspace.  In  this  paper  similar  techniques  are  used  to 

j obtain  an  analytical  solution  to  the  problem  of  the  propagation  by  a 

rigid  wedge  of  a 2-dlmenslonal  crack  in  an  Infinite  power  law  material. 

This  particular  problem  Is  an  example  of  a physically  realizable 
mode  of  fracture  to  which  the  correspondence  principle  does  not  apply. 
Moreover,  the  closed  form  solutions  obtained  here  provide  a method  for 

■ experimentally  determining  physical  parameters,  such  as  the  stress 

•» 

intensity  factor.  Important  for  an  understanding  of  fracture  phenomena. 

This  problem  has  been  studied  extensively  for  elastic  material,  a 
detailed  description  of  which  appears  in  the  article  by  Barenblatt  [2]. 
The  wedging  problem  was  considered  by  Barenblatt  for  two  reasons.  In 
general,  Barenblatt  addressed  himself  to  the  Issues  of  unrealistic 
singular  crack  tip  stresses  and  Indetermlnant  crack  lengths  encountered 
In  fracture  In  classical  linear  elasticity.  The  wedge  problem  exhibits 
both  features  even  though  the  crack  Is  "semi-lnf Inlte"  since  the  free 
crack  surface  ahead  of  the  wedge  Is  finite.  In  particular,  the  thrust 
of  Barenblatt 's  article  Is  that  In  a neighborhood  of  the  crack  tip, 
the  separated  surfaces  experience  cohesive  molecular  forces  which 
produce  a cusp  profile  and  a finite  stress  field.  Barenblatt  accounts 
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for  this  by  the  phenomenological  expedient  of  assuming  that  for 
equilibrium  cracks  a cohesive  normal  traction  acts  over  a small 
Interval  behind  the  crack  tip  with  a constant  stress  Intensity  factor 
characteristic  of  the  material.  As  noted  by  Barenblatt,  the  finite 
crack  length  is  determined  by  merely  knowing  this  stress  Intensity 
factor  and  Is  independent  of  the  actual  distribution  of  cohesive  forces. 

The  utility  of  this  observation  lies  In  reversing  the  implication. 

Namely,  if  experimentally  a wedge  Is  driven  Into  an  elastic  material 
and  the  length  of  the  free  surface  of  the  resultant  equilibrium  crack 
is  measured,  then  the  stress  intensity  factor  may  be  easily  calculated. 

Here  we  provided  simple  closed  form  expressions  for  the  stress 
and  displacement  fields  for  the  wedging  of  a power  law  viscoelastic 
material  under  the  Barenblatt  hypothesis.  This  analysis  culminates 
In  a relation  among  the  material  parameters,  the  material  stress 
Intensity  factor  and  the  length  of  the  free  crack  surface. 

Mathematically,  the  wedging  problem  reduces  to  a set  of  dual 
Integral  equations  which  In  a sense  are  Inverse  to  those  considered  In 
WNS.  These  equations  are  first  solved  without  the  Barenblatt  hypothesis 
to  simply  see  what  the  stress  and  displacement  fields  are  for  a power 
law  viscoelastic  material. 

It  was  observed  that  the  nature  of  the  crack  tip  profile  is 

■ 


dependent  on  the  material  whereas  the  nature  of  the  stress  singularity 
Is  not.  In  particular,  the  classical  square  root  singularity  In  the 


normal  stress  occurs  for  all  materials,  whereas  crack  tip  profiles 
range  from  near  square  root  (as  in  the  elastic  case)  to  cusp-like. 
Hence,  singular  stresses  can  occur  even  in  the  presence  of  a smoothly 
closed  crack.  With  the  addition  of  the  Barenblatt  hypothesis,  the 
unknown  crack  length  is  then  determined. 


SECTION  1.  FORMULATION  OF  THE  PROBLEM. 

We  now  consider  the  problem  of  the  steady  translation  to  the  left 
(with  velocity  U)  of  a rigid  2-dlmenslonal  S5nnnnetrlc  wedge  with  finite 
asymptotic  thickness  In  an  Infinite  power  law  viscoelastic  body.  (See 
Fig.  1.)  Due  to  S3mimetry,  we  may  consider  the  problem  for  the  upper 
half  plane  only.  The  notation  adopted  here  Is  consistent  with  that 
employed  In  WNS  [1]. 

Neglecting  Inertia,  the  force  balance  equations  are 

9o  3t 

9x  3y  ^ 


3t 

+ ^ 


and  the  boundary  conditions  are 


Txy(x,  0,  t)  = 0 


v(x,  0,  t)  - 


i " 

(f(x- 


-00  < X < ® 


X < Ut 


a^(x,  0,  t)  - 0 


Ut)  X > JI2  + Ut 


Ut  < X < ^2  + Ut. 


Here  v(x,  y,  t)  and  u(x,  y,  t)  are  the  vertical  and  horizontal 

displacements,  a and  a are  the  normal  stresses  on  x^constant 
X y 
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lines  and  y=constant  lines  respectively,  ^ is  the  shear  stress 

and  f(x)  is  the  wedge  profile  whose  finite  asymptotic  thickness  is  h, 

i.e.  lim  f(x)  = h.  The  viscoelastic  stress  strain  laws  are 
x-»  “ 


■ f 
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A(t-T)  di  + 2 


A(t-T)  1^  dx  + 2 


ft  9e 

G(t-T)  ^ dx 


ft 


9e 
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G(t-,)  ^ d. 
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the  strains. 

A detailed  description  of  the  next  series  of  steps  appears  in  WNS 
[1]  and  will  be  omitted  here.  After  adopting  the  Gallilean  variable 
s = X + Ut,  all  of  the  equations  are  differentiated  with  respect  to  s 
(save  for  the  C7y(s,  0)  boundary  condition)  and  then  Fburier  transformed 
with  respect  to  s.  The  resulting  ODE's  in  y are  solved  and  a Fourier 
Inversion  results  in  the  following  integral  equations  relating  v'(s,  0) 
to  the  normal  traction  Oy(s,  0) 


-r 

^ -a 

-r 


. .1 


e^^’^Aj  (p)dp 


dp 


Oy(s,  0) 


v’(s,  0) 


-»  < s < 


.00  < 3 < 00 


where 


3o 

■57^  (s,  0) 


2(27r)^^^ip)|p|  A (p). 

(A+2G) 


(Here,  F denotes  the  Fourier  transform  of  F. ) 


To  characterize  our  power  law  material  It  suffices  to  assume  that 

G(A-fG)  ^ E 
(A+2G)  4(l-v^) 

where  E(t)  = Et  ^H(t)  and  v Is  Poisson's  ratio  which  Is  taken  to 
be  constant.  This  leads  finally  to  the  dual  relations 


1 


dp 


v'(s,  0) 


►00 

Aj^(p)  |p1  (Ip) 

-*00 


1 Isp 
e 


dp 


g(s) 


-00  < g < 00 


where 

g(s)  = 2(l-v^)0y(s,  0)/(u“Er(l-a)). 

It  should  be  noted  that  here  v'(s,  0)  Is  known  off  of  the  compact 
Interval  (0,  on  which  g(s)  Is  known.  Hence,  the  dual  relations 

(1)  and  (2)  are  Inverse  to  those  derived  In  WNS  [1]. 

The  first  problem  considered  Is  that  of  finding 
v'(s,  0)  (=  v'(s)  for  0 < s < and  where  0 < < ^2 
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and  = 0,  given  that  g(s)  = 0 for  0 < s < 42’  As  will  be 

observed  later,  this  problem  is  underposed  in  that  (the  length  of 

the  free  crack  surface)  and  (the  distance  from  crack  tip  to  wedge 

point)  cannot  both  be  determined.  Indeed,  only  an  explicit  value  for 
the  ratio  emerge. 

To  solve  (1)  and  (2)  we  first  rewrite  (2)  as 

r Aj^(p)  IpI  (ip)“  ^e^^^dp  = g(s). 

/ZlT 


The  representation  of  Aj^(p)  obtained  from  (1)  suggests  the  validity 
of  the  formal  manipulations  leading  to  line  (3).  The  inversion  of  (1) 
and  substitution  of  the  resulting  expression  for  Aj^(p)  into  (3) 
yields  the  following  integral  equation 


/CO 

d 

V 

ds 

J .00 


(y)K(s-y)dy  = g(s) 


-00  < s < 00 


where 


K(z)  = - 


i(a7r)  |z| 


z > 0 


z < 0. 


The  operator  appearing  on  the  left  side  of  (4)  is,  in  a generalized 
sense,  inverse  to  the  Abel  type  operator  discussed  in  WNS  [1].  This 
point  requires  some  ampllficntlon.  If  denotes  the  operator  of  WNS 

[11  on  say  (0,  , i.e. 


/ 


K^(f) 


Kj^(s-y)f  (y)dy 


K^(z) 


r(i-a) 

TT 


-cos(a'iT)  I z 
la-l 


la-l 


i-r 


z > 0 
z > 0, 


and  if  K denotes  the  operator  in  (4)  then  Samko  [3]  proved  that  the 
product  KK^  is  a weighted  finite  Hilbert  transform  which  in  a certain 
sense  acts  like  an  identity  operator.  On  (-«>,  “)  and  K are 
true  inverses.  Hence,  if  either  f'  or  g is  known  on  all  of 
(_oo^  oo)  j the  other  may  be  easily  determined.  However,  since  here  and 
in  WNS  [1]  mixed  boundary  conditions  are  assumed,  no  simple  Inversion 
of  the  generalized  Abel  equations  exists. 

At  this  point  it  is  necessary  to  observe  that  the  nature  of  the 
wedging  problem  changes  dramatically  for  o > 1/2.  Indeed,  from  the 
analysis  of  the  asperity  problem  contained  in  WNS  [1],  it  is  apparent 
that  the  material  will  not  remain  in  contact  with  a wedge  of  finite 
asymptotic  thickness,  but  rather  will  leave  off  at  some  point  with 

^3  ^ ^2‘  Fig"  !•)  Consequently,  the  appropriate  boundary 

conditions  for  a > 1/2  are 


v(s,  0) 


0 s < 0 

f(s)  ^2  s < 


I 


3 


< 


s. 


Oy(S,  0)  = 0 


0 < s < 
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In  this  case,  solving  equation  (4)  is  just  as  difficult  as  solving  the 
multiple  asperity  problem  mentioned  in  WNS  [1],  (Section  3).  If  the 
wedge  profile  is  monotone  increasing,  this  pathology  does  not  occur  for 
0 ^ a _<  1/2.  Therefore,  detailed  solutions  will  be  presented  for  this 
case  only.  Even  though  an  explicit  solution  for  a > 1/2  cannot  be 
exhibited,  it  will  prove  to  be  an  easy  matter  to  infer  from  the 
solution  for  0 < a < 1/2  what  the  nature  of  the  stress  distribution 
and  crack  profile  is  in  the  vicinity  of  the  crack  tip  when  a > 1/2. 

Henceforth,  we  shall  assume  that  0 < a < 1/2,  since  the 
limiting  cases  a = 0 (the  classical  elastic  case)  and  a = 1/2  may 
be  then  easily  deduced.  The  case  a = 1/2  may  also  be  treated 
directly  in  a simple  manner,  since  the  integral  equation  (4)  reduces 
to  an  ordinary  Abel  equation. 

To  solve  (4)  when  0 < a < 1/2,  we  first  integrate  and  recall 
that  on  (-“>,  0),  v' (s)  vanishes  and  on  (.^2^  “)  > v' (s)  = f'(s),  a 
known  function.  Therefore  we  obtain  the  generalized  Abel  equation 


where 


r(a) 

IT 


'^2  v’ (t)dt 

/ X Ct 

_s  (t-s) 


+ 


cos(aiT) 


s 

0 


v' (t)dt 
(s-t)“ 


F(s) 


0 < s < 


F(s) 


r(a)  r v’(t)dt 


(5) 
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and  C is  an  arbitrary  constant  of  integration.  Equation  (5)  is  the 
equation  considered  in  WNS  [1]  with  oc  replaced  by  1 - a.  Its 
solution  may  be  obtained  either  by  the  methods  employed  there  or  by 
that  of  Samko  [3].  We  adopt  the  former  approach.  However,  the  details 
of  the  derivation  of  the  solution  will  be  omitted  since  the  essential 
ideas  are  contained  in  WNS  [1]. 

We  remark  that  for  technical  reasons  a straight  substitution  into 
the  formulas  derived  in  WNS  [1]  cannot  be  used  to  solve  the  example 
considered  in  the  next  section  since  it  results  in  divergent  integrals. 

It  proves  more  convenient  to  solve  the  corresponding  Rlemann-Hilbert 
boundary  value  problem  in  the  class  of  functions  unbounded  at  both 
s = 0 and  s = 8.2,  thereby  introducing  an  extra  arbitrary  parameter, 

K,  which  is  later  specified  by  imposing  continuity  conditions.  The 
details  will  be  omitted  and  only  the  final  results  presented. 

It  may  be  shown  by  the  method  employed  in  WNS  [1]  that  v'(s) 
must  satisfy  the  ordinary  Abel  equation 

= H.  (s ) 0 < s < 8 (6) 

•'0  (s-t)  ^ 


with 


H^(s)  = Ks"^'^^(82-s)^'^^““ 


-1/2 


(7) 


i 

f 

! 


The  integral  in  (7)  is  a Cauchy  principal  value. 
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ft 


The  function  v(s)  is  easily  determined  from  (6)  to  be 


v(s) 


^ Hj^(t)dt 

• Jo  (s-t)^-“ 


B(a 


(8) 


where  B(*,  •)  denotes  the  beta  function.  The  parameters  K,  C and 
^2^^1  then  chosen  so  as  to  insure  the  continuity  of  v(s)  at 

s = 0,  and  v'(s)  at  £2* 

In  the  next  section,  an  example  is  considered  for  which  a simple 
closed  form  solution  is  obtained  for  v(s).  We  close  this  section  with 
a useful  observation.  Once  v(s)  is  determined,  g(s)  for  s < 0 
is  easily  calculated  from  (4).  In  particular,  we  have 


g(s) 


IM 

IT  ds 


»00 

0 


v'(v) 

(y-s)“ 


dy. 


(9) 


It  may  be  shown  from  (7)  that  K and  C can  be  chosen  so  that 

v(s)  = + v^(s)  (10) 

where  Vj^'(s)  is  continuous  at  s = 0.  Hence  for  0 ^ a < 1/2, 

v(s)  has  a vertical  tangent  at  the  crack  tip,  s = 0,  and  from  (9) 

—1/2 

it  follows  that  g(s)  = 0((s|  ) as  s approaches  zero  from  the  left. 

It  is  apparent  that  (10)  is  also  valid  for  1/2  £ a < 1. 

Consequently,  for  a = 1/2  the  free  crack  faces  meet  at  a sharp  corner 
and  there  results  a square  rout  singularity  in  the  normal  stress  ahead 


i 


i 


of  the  crack  tip.  However,  for  a > 1/2,  the  crack  surfaces  form  a 
cusp  (i.e.  close  smoothly)  with  the  normal  stress  still  having  a 
square  root  singularity.  (See  Fig.  1.)  This  is  in  marked  contrast 


to  the  behavior  of  elastic  materials. 
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SECTION  2.  AN  EXAMPLE. 

To  illustrate  the  technique  of  the  previous  section  we  choose  a 
wedge  shape,  f(s),  for  which  the  Cauchy  integral  in  (7)  is  easily 
computed.  Specifically,  for  0 £ a < 1/2  we  take 

f(s)  = h[l  - (s/J,^)®  ^] , s The  fact  that  the  wedge  profile 

depends  on  a does  not  impose  a serious  limitation  on  the  usefulness 
of  the  example,  since  the  shape  of  f(s)  is  qualitatively  the  same  for 
all  values  of  a. 

With  this  choice  of  f(s),  the  function  F(t)  in  (5)  becomes 
F(t)  = C - Jl^^"“hT"^[l  - (I-t/A^)^"®]. 


The  Cauchy  integral  in  (7)  may  now  be  evaluated,  and  after  some 


manipulation  yields 


0 


F(T)dT 


a-1/2 


= C B(a+l/2,  l/2)i2  (l-s/)l2)  2^1^^^^'^'  1/2;  s/i^) 


Xh  !l(.“tl/2j.r(-l/2)  (i_s/ii^)“-1/2^F^(1/24<x,  1/2;  3/2;  s/l^) 


- Xhr(a), 


where  the  parameter  A is  defined  by 
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X = (12) 


After  substitution  of  (11)  into  (7),  the  integral  in  (8)  may  be 
computed  from  which  it  follows  that 

v(s)  = (K+Xh)  1/2;  o+l/2;  s/l^) 


— — (s/Jl  - s/l  )^^^ 

r(a+l)  '^2' 


\h2{s/l^f~^^^^?^{l/2,  1/2;  3/2;  S/JI2). 


For  convenience,  we  define  a new  variable  z by  z * 3/^2 
consider  v(s)  as  c function  of  z,  denoted  v(z).  To  determine  C, 
K and  X we  Impose  continuity  for  v(z)  near  z = 0,  1 and  for 
v'(z)  near  z = 1.  It  is  readily  seen  that  near  z = 0 v(z)  has 
the  following  form 


A,  , „ a-l/2  a+l/2  a+3/2v 

v(z)  = C^z  + C2Z  + 0(z  ), 


(13) 


where 


CitA^ 


C,  = B(o,  1/2)  (K+Xh)  - 


r(a+l) 


and 


- 2Xh 


Cit4, 


(14) 


C2  - B(a,  1/2)  2(a+l/2)  2r(a+l)  " 2^' 


(15) 
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A similar  analysis  near  z » 1 yields 

v(z)  = + D2(1-z)^^^  + 0((l-z)^^^) 


(16) 


with 


D 


1 


» K 


(17) 


and 

- (K+Xh)B(l/2,  a)(l-2a)  - + 2Xh.  (18) 


Continuity  of  v(z)  at  z = 0 and  z «■  1 require  that  Cj^  - 0 and 
^1  * ^^^2^  “ h(l-X),  and  continuity  of  v' (z)  at  z = 1 follows  from 
D2  = 0.  Hence,  we  obtain  the  relations 


nCi 

2_  _ r(a+l)r(l/2)h(l-ot) 

r(a)  " r(a+l/2) 

, . r(3/2)r(a-H) 

r(a+l/2)  • 

The  constants  K and  X are  determined  explicitly,  whereas  C and 
I2  occur  coupled  In  relation  (19).  The  addition  of  the  Barenblatt 
assumption,  discussed  In  the  next  section,  will  provide  the  required 
extra  relation  for  the  determination  of  C and  l2* 

It  should  be  noted  that  the  coefficient  C2  becomes 


(19) 


(20) 
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C 


2 


2r(l/2)r(a+l)(l-a)h 
r(3/2+o)  3 


4 

3 


Xh 


■ a-<>) 

(1/2+a) 


(21) 


vhlch  is  clearly  non-zero.  This  verifies  line  (10)  In  section  1 and 
hence  the  subsequent  discussion  on  the  nature  of  the  singular  nomal 
stress. 


I 

I 
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SECTION  3.  THE  BARENBLATT  HYPOTHESIS. 

The  Barenblatt  hypothesis  Is  that,  for  an  equilibrium  crack, 
forces  of  cohesion  act  on  the  separated  crack  surfaces  within  a small 
neighborhood  of  the  crack  tip  to  effect  a smoothly  closed  crack  and  a 
cancellation  of  the  singular  normal  stresses  ahead  of  the  crack  tip. 

It  Is  further  assumed  that  the  stress  Intensity  factor  corresponding 
to  these  cohesive  forces  Is  an  intrinsic  property  of  the  material 
and,  in  particular.  Is  Independent  of  the  length  of  the  Interval  on 
which  the  forces  act.  The  Barenblatt  hypothesis  Is  encorporated  into 
the  model  discussed  In  section  2 by  amending  the  stress  free  normal 
boundary  condition  as  follows: 

0)  “ p(s)  0 < s < 

where  p(s)  Is  some  distribution  of  cohesive  normal  tractions  which 
vanishes  for  s d and  d is  a small  positive  number.  The  tractions 
p(s)  are  assumed  to  have  a known  stress  intensity  factor,  L,  given 
by 

I,  . 1 r*  a,. 

' Jo  ^ 

It  Is  clear  that  the  dual  Integral  equations  (1)  and  (2)  are  still 
valid.  However,  the  generalized  Abel  equation  (5)  becomes 
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where 


r(a) 

n 


v' (t)dt 
(t-s)“ 


fS 


+ cos (an) 


0 


v'(t)dt 


= F^(s) 

0 < s < 


F.  (s)  = F(s)  - Q(s), 


Q(s) 


q(t)dt 


0 < s < d 


and 


q(s)  - 2(l-v^)p(s)/(u“Er(l-a)). 


Equation  (22)  may  be  reduced,  as  before,  to  the  ordinary  Abel 
equation 


r - H2(s). 

Jo  (s-t)®  ^ 2 


where  Hj^(s)  Is  given  by  (7)  and 
H (s)  - -s"J^^^(!l  -s) 


Solving  (26)  we  obtain,  finally. 


v(s) 


B(a,  1-a) 


.1/2-a  ril-_a)  f^2 

’ Jo 

(t-s) 

Hj(t)dt 

fS  H2(t)dt 

jo  (s-t)^““ 

JQ  (8-t)^"“ 

Q(T)dT 


(22) 

(23) 

(24) 

(25) 

(26) 

. (27) 


(28) 
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A few  comments  are  now  In  order.  As  Is  well-known  (Barenblatt  [2]), 
it  Is  Che  stress  Intensity  factor,  L,  rather  than  the  actual 
distribution  of  tractions,  p(s),  that  enters  Into  an  analysis  of 
stresses  at  the  crack  tip  and  a determination  of  the  free  crack  length, 
Ij-  Moreover,  we  shall  show  that.  If  is  small,  then  the 

displacement  field  is  essentially  independent  of  p(8).  Consequently, 
for  Illustrative  purposes,  we  shall  assume  that  p(s)  Is  equal  to  a 
constant,  p,  on  0 < s < d. 

Returning  to  Che  computation  of  v(s),  we  denote  the  constant 
function  q(s)  In  (25)  by  q and  observe  that  Q(s)  In  (24)  becomes 


Q(s)  - q(8-d) 


0 < s < d. 


Substitution  of  Q(s)  Into  (27)  yields 


H r N - .1/2- 

H2(s)  “ - s (i2"s) 


a r 


d T^/^(i2-T)“‘^^^(T-d) 


(t-s) 


dr. 


O ^ / o 

Due  to  the  presence  of  (12"’^)  the  Integrand,  the  Cauchy 

Integral  In  (29)  cannot  be  evaluated  simply.  Since  we 


(29) 


t 
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ultimately  will  let  AII2  tend  to  zero,  we  may  substitute 
for  in  (29). 

Given  this  simplification,  H2(s)  becomes 

/ X / /jn-1/2,,  /„  ^l/2-o  r(l-a) 

H (s)  = (s/d)  (i  - s/Jl„)  q — ^ 

1 Z 7T 

d B(2,  1/2)^?^ (-3/2,  1;  1/2;  s/d)  0 < s < d 
-d(s/d)~^B(2,  3/2)2F^(1,  3/2;  7/2;  d/s)  d < s. 

To  analyse  v(s)  in  (28)  we  introduce  some  notation.  Let  Vj^(s) 
denote  the  displacement  corresponding  to  a free  crack  surface,  i.e. 
without  the  Barenblatt  hypothesis.  Then 


Vj^(s)  - 


^ j-s  H^(t)dt 
'0  (s-t)^-“ 


B( 


0 < s < £, 


and  is  given  explicitly  by  (13).  Let  V2(8)  be  given  by 


V2(s) 


1 


/■s  H2(t)dt 


B(a,  1-a) 


0 (s-t) 


1-a  • 


As  in  the  preceedlng  section,  it  is  convenient  here  to  define  new 
variables  z - and  e = and  let  v,  Vj^  and  V2  denote  the 

corresponding  functions  of  z and  e.  Moreover,  we  Introduce  the 


parameter 
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L 


1 


L 2(l-v^)/(v“Er(l-a)) 


where  L is  the  stress  Intensity  factor  corresponding  to  the  cohesive 
tractions  p(s)  on  (0,  d). 

It  is  now  a simple  matter,  albeit  tedious,  to  substitute  the 
expression  for  H2(s)  into  (28)  and  compute  the  necessary  Integrals, 
keeping  in  mind  that  H^Cs)  is  defined  differently  on  (0,  d)  and 
(d,  ' However,  our  ultimate  goal  is  the  determination  of  and 

^2  for  which  we  need  only  consider  the  asymptotic  form  of  ^(z)  for 
z near  0 and  1.  Moreover,  we  are  assuming  e to  be  small,  and 
hence  shall  combine  all  terms  which  vanish  as  e -»•  0 Into  a single 
residual  term. 

We  first  consider  the  asymptotics  of  v(z)  for  z near  zero,  in 
particular  we  take  z < e.  It  may  be  shown  that 


V2(z)  = (2L^/3r(o))[eB(o,  1/2) £2'^^^^z“  1/2;  a+1/2;  z) 

- 3B(a,  3/2))l2'^^''^z“'^^''^2^l(“-l/2.  3/2;  3/2-hi;  z) 


^ l/2„ 

+ E £„ 


1/2^  r 

Jo  (z-t)^"“  ^ ^ 


0 < z < e 


(30) 
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with 


and 


C3  = (2L^/3)e£2'*’^^^r(l/2)/r(o+l/2) 


= l/2C3(a-l/2)/(a+l/2)  - Lj^z“'^^^^r(l/2) /r(a+3/2) . 


Since  v(z)  = Vj^(z)  - we  obtain  for  0 < z < e 

v(z)  * (C^-C3)z“"^''^  + + 0(z“''’^''^). 


For  finite  stress  we  require 


C3  - C3  = 0 


C2  - = 0. 


A similar  analysis  must  be  performed  for  z near  1.  Observe 
first  that,  if  e < z < 1, 


>>2(2)  * V3(z)  + v^(z) 


with 


'2  r 

' ■ BCc  1-a)  J„ 


and 

“ B(a,  1-a)  ‘ 

Clearly  v^Cz)  is  smooth  for  z near  1 and  Is  negligible  for 
e « 1.  For  v^(z)  we  obtain 

^^(z)  - D3  + D^(l-z)^^^  + 0((l-z)^/^) 

with 

D3  = - 0(e).  (34) 


3/2 

Since  (e/t)  3/2;  7/2;  e/t)  is  uniformly  bounded  for 

e £ t ^ z _<  1,  the  estimate  (34)  is  now  evident. 

For  continuity  of  v(z)  and  v'(z)  at  z * 1 we  require 

- D3  - f(Jl2>  “ 

and 


(35) 


°2  ■ °4  “ 
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Combining  (31),  (32),  (33),  (34)  and  (35),  we  arrive,  finally,  at  the 
relations 

K = h(l-X)  + 0(e), 

X = r(3/2)r(l+a)/r(l/2+a)  + 0(e), 

CTDl® 

-j~  = r(l+a)r(l/2)h(l-a)/r(l/2+a)  + 0(e) 

and 

^a+1/2  ^ _(2/3)h(l-a)r(]+t.)/Lj^.  (36) 

We  remark  again  that  line  (36)  provides  a means  of  experimentally 
determining  the  stress  Intensity  for  an  equilibrium  crack.  Indeed, 
from  (36)  It  follows  that  the  stress  Intensity  factor,  L,  associated 
with  the  cohesive  tractions  (which,  under  the  Barenblatt  assumption. 

Is  equal  In  magnitude  but  opposite  In  sign  to  the  usual  stress  intensity 
factor  for  a free  crack  surface)  is  related  to  the  wedge  speed,  U, 
and  crack  length,  by 

L = -(l/3)hr(l+a)r(2+a)u“Ef~““^^^/(l-v^). 

Hence,  for  given  speed  U,  the  stress  Intensity  factor  L Is  easily 
computed  once  the  length  of  the  crack,  I2*  has  been  measured. 


ji 
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A few  comments  on  the  stress  and  displacement  field  under  the 
Barenblatt  hypothesis  are  In  order.  Outside  of  a neighborhood  of 
order  e centered  at  the  crack  tip,  the  stress  and  displacement  fields 
are  merely  a perturbation  of  order  e from  the  stress  and  displacement 
fields  corresponding  to  a free  crack  surface  (l.e.  without  a 
Barenblatt  zone) . It  is  only  on  an  Interval  of  order  e that  the 
displacement  profile  suddenly  becomes  cusp-llke  and  the  stress 
distribution  falls  off  rapidly  to  zero.  In  the  limit  e = 0,  (though 
keeping  L constant)  the  stress  and  displacement  fields  are  Identical 
with  those  of  a free  crack,  every  where.  However,  unlike  the  problem 
posed  with  a free  crack  surface,  we  are  still  able  to  determine  the 
crack  length,  explicitly. 
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